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Orbital period derivative of a binary system using an exact orbital energy equation 
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It is proposed that the equations of motion which yielded two major predictions of general relativ- 
ity without utilizing Riemannian geometry and geodesic trajectories are exact in nature and can be 
applied to pulsars and inspiralling compact binaries for analyzing orbital period derivative and two 
polarization gravitational wave forms. Exactness of these equations eliminates the need for higher 
order iPN corrections to the orbital energy part of the balance equation. This is mainly due to the 
introduction of dynamic WEP which states that the gravitational mass is equal to the relativisitc 
mass. 
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1. INTRODUCTION 

The famous quadrupole formalism of Einstein Q is the 
lowest order wave generation formalism in the Newtonian 
limit 1/c — > 0. Using this formalism, Peters and Mathews 
[Tol | calculated the time average of the energy radiated 
from a system of two gravitating masses. From these re- 
sults, secular decays of the semimajor axis, and eccentric- 
ity were found by Peters Q • Based on this work, Esposito 
and Harrison [ll[ , and Wagoner [l^] heuristically formu- 
lated the orbital period derivative for the binary pulsar 
system in light of the discovery of a radio pulsar PSR 
1913+16 by Hulse and Taylor pj|]. Timing measurement 
of PSR 1913+16 by Weisberg and Taylor [3 provided 
confirmation for the rate of orbital period decay due to 
gravitational radiation damping as predicted by general 
relativity. Subsequently both the theoretical model and 
the measurement accuracies were im pro ved over a period 

of time H EE EE EE EE EE HE HEM M ■ 

In periodic relativity theory (PR) [l| , we extended the 
concept of space-time symmetry to include heavier de 
Broglie particles, which enabled us to introduce energy- 
momentum invariant into the space-time invariant result- 
ing in an invariant relationship between the force and 
the energy. We utilized this relation to address two-body 
problem in gravitational field by proposing the following 
equation of orbital motion 



dt 



(E - m c 2 ) 



um^dr , . . 
— 2-r— (cosip + smtp) , (1.1) 
r z at 



this set-up for the equation of motion, we successfully 
obtained exactly the same lowest order expressions for 
gravitational redshift and deflection of light as in gen- 
eral relativity. We also obtained second order non- 
homogeneous non-linear differential equation of motion 
equivalent to general relativity which yielded same val- 
ues of the pcrihclic precession for the planets of the so- 
lar system. These are reasons enough to conclude that 
Eq. (jl.ip is an exact equation of motion Q and can pro- 
vide us exact orbital energy for two body systems. If 
we indiscriminately introduce a;PN corrections to these 
equations, there is a good chance that we may not be 
able to derive correct formula for gravitational redshift 
and bending of light from such corrected equations. This 
means that we have locked on to an orbital energy equa- 
tion which we cannot alter. This is mainly because the 
present theory is energy based theory. In comparision, 
general relativity uses Newtonian orbital energy equa- 
tions and adds hi ghe r ord er p ost-Newtonian corrections 
[E H HI HE 0, 12E EE S_f to get desired relativistic 
effects for computing center of mass binding energy and 
gravitational wave energy flux. In PR, these corrections 
are fundamentally built into the orbital energy equation. 
The basis of this equation is the introduction of dynamic 
WEP (weak equivalence principle) in PR. Dynamic WEP 
states that the gravitational mass is equal to the rela- 
tivistic mass. For more details see [l|. This situation 
provides an ideal testing ground for proving this theory 
in the area of gravitational radiation. 



where m — 77710 is the relativistic mass, E total energy 
of the particle, if> is the angle between the radial vec- 
tor r and the velocity vector v. Factor (cosip + sxnip) is 
responsible for giving us geodesic like trajectories. Ap- 
pearance of this factor is due to the introduction of the 
Newtonian gavitational potential in the form of summa- 
tion of two components due to acceleration (/xm/r ) cos ip 
and (/im/r 2 ) sin ip acting on the gravitating body. With 



ORBITAL ENERGY OF A BINARY SYSTEM 



For a binary system, we can rewrite Eq. (jl.ip as 



^-(E-m 2 c 2 ) = - Ml "! 27 v(cos7/; + sinV')f, (2.1) 
dt r z 



'Electronic address: cons_engl@yahoo.com 



where mi is the mass of the pulsar and 777,2 that of the 
companion star. For Binary, /ii = Gm\, r is the separa- 
tion. Total orbital energy E of the system can be given 



2 



by 



E° = (E - m 2 c 2 ) + J 



/iiTO 2 7 dr 
r 2 dt 



(cos ip + s'mtp) dt. 



(2.2) 



If we evaluate E at periastron, we can intorduce fol- 
lowing approximations, ip = tt/2 and 77127 = constant. 
Therefore, 

E° = (E-rn 2 c 2 ) + W m 27 (-- ) , (2.3) 



E° = m 2 c 2 



(7 - 1) 2" 

rcr 



£° = m 2 c 2 



c 2 ) i rc 2 } 



(2.4) 
(2.5) 



where v p is the velocity at periastron. One can observe 
that the xPN corrections are naturally and fundamen- 
tally built into the orbital energy equation Eq. (12.5[) . 
Unlike general relativity, there is no limit to which this 
equation can be expanded. However, here we will expand 



Eq. (|2.5p equivalent to 3.5PN order in general relativity. 



E° = m 2 c 2 



1 



13. 

2 c 2 



3«2 



5 vl 35 vl 
16 c6 + 128 



E u 



8 c 4 

{ rc 2 1 

Him 2 1 /J 2 m 2 (l + e)(l - 3e) 
8 c 2 a 2 (l -e) 2 

1 /i?m 2 (l + e) 2 (l - 5e) 
+ 16 c 4 a 3 (l-e) 3 

5 n\m 2 {l + e) 3 (l - 7e) 



128 



c 6 a 4 (l -e) 4 



(2.6) 



(2.7) 



where v p = /i/r p , ft 2 = (i%a(l — e 2 ) and r p = a(l — e) 
at periastron. Variation of £7 with time is then given 
by Eq. (|2.8j) . It is to be noted here that introduction of 
the deviation to flat Minkowski metric [J in Eq. |(23J) 
and Eq. (|2.8p does not alter the formulation because the 
effect gets cancelled out. 



dE° 
dt 



y>\m<x 1 /j 2 to 2 (1 + e)(3e - 1) 
2a 2 ~~ I" 



3 /x?m 2 (l + e) 2 (5e - 1) 5 u\m 2 {\ + e) 3 (7e - 1)" 



/i 2 m 2 e 
c 2 a 2 (l — e) 3 



c 2 a 3 (l — e) 2 
1 u\m 2 (7e 2 



4e 



16 
1) 



c 4 a 4 (l - e) 3 
15 fx\m 2 e{l + e) 2 



32 



c 4 a 3 (l - e) 4 



C 6 a 4(!_ e )E 



O 



c 6 a 5 (l 
1 



e) 4 



da 
lit 



(2.8) 



ORBITAL PERIOD DERIVATIVE 



nology by 



Laser interferometric observations of gravitational wa- 
ves demands energy balance equation computed to an 
extraordinarily high degree of accuracy of order 1/c 6 . 
Equations of variation with time of the orbital frequency 
and orbital phase of an inspiralling compact binary are 
derived from this energy balance equation and the the- 
oretical templates of the compact inspiral binary is ob- 
tained by introducing these highly accurate parameters 
into binary's two polarization wave-forms h + and h x 
In case of binary pulsars, the energy balance equation is 
given by 



dE° 
dt 



= -£, 



(3.1) 



where E is the orbital energy and £ the total gravita- 
tional luminosity (or wave flux) of the source. The time 
average value (£) has been computed to 1PN order by 
Blanchet and SchaferQ and can be given in our termi- 



_ 1024 v 2 {-E»f , , 

W ~5GcV(l- 
-E° 



73 



37 



— -I e* 

24 96 



+ 



168/_tc 2 (l -e 2 ) 
6419 



-840 



13 - 6414^ 
, 5103 . 



27405 



5377 



259 



8 



e 6 I v 



(3.2) 



where /1 = (niim 2 )/m, m = (mi + m 2 ), v = /i/m and 
E° to 1PN accuracy can be obtained from Eq. (|2.7|) and 
given by 



V\m 2 
2a 



1 | l m (l + e)(l-3e) 
4 ac 2 (l-e) 2 



(3.3) 



We can truncate Eq. (|2.8|1 to 1PN accuracy as follows. 
dE° 



dt 



uim 2 _ 1 /z 2 m 2 (l + e)(3e - 1) 
2a 2 ~ 4 



/-i 2 m 2 e 



c 2 a 2 (l — e) 3 



c 2 a 3 (l — e) 2 
de 
dt' 



da 
~dt 



(3.4) 
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We will utilize following lowest order expression for 
(de/dt) given by Peters Q. 1PN expression could be 
more consistent but it may not affect the result signifi- 
cantly. 



de 



19e(l - e 2 )(f + (f2f/304)e 2 ) 



dt I 12a(l+ (73/24)e 2 + (37/96)e 4 
We get following relation from Kepler's third law, 

da 2a Pb 
dt = ift' 



(3.6) 



where Pb is the orbital period and Pb orbital period 
derivative. Substitution of Eqs. (|3.5p and (|3.6p in 
Eq. 1331) yields, 



\ dt I ~ 3a 



Hi(ai + a 2 ) 



CT3 



Pb 
P b 



(3.7) 



where cri = (-3 + 6e + (151/8)e 2 + (149/4)e 3 ), 
ct 2 = ((1081/32)e 4 + (79/8)e 5 + (lll/32)e 6 ), 
cr 3 = 6ac 2 (l - e) 2 (l + (73/24)e 2 + (37/96)e 4 ). 



Introducing Eqs. (|3.2|) and (|3.7p into the balance Eq. (|3.1 
and replacing a with Kepler's third law in the form 



f 

a 



we get, 



/i 2 3 , where fi 2 = G(mi + m 2 ), 

(3.8) 

ft = fci(l + fe 2 ), where, (3.9) 



1927T / 2nG\ 3 m\m 2 (m\ + m 2 ) a 
(l-e 2 )i 



5c 5 V p b 

■ , 73 2 37 4 

x ^ H e 2 H e 4 

' 24 96 



fc 2 = 



fafafafa 



(3.10) 



(3.11) 



47T 
= - — 



rnim 2 (mi + m 2 ) 3 /2ttG\' c 
"ft 



fc4 = 



'!•") c 7 (l - e 2 ) 
1 (l + e)(l-3e) {2ttGY 



(3.13) 



1 + 4 c 2 (l-e) 2 



ft, 



mi(mi + 3 
(3.14) 



fa = 



13 -6414e*-^e 4 -^ 
4 16 



-840 



6419 2 5103 4 
p z p. 4 - 



259 6 \ 
J" 



(3.15) 



For binary pulsar PSR 1913+16, factor k 2 given by 
Eq. (|3.1ip in this theory turns out to be k 2 = —1.865 x 
10~ 5 compared to +2.15 x 10 -5 given by Blanchet and 
Schafer 0], which is —1.1528 times factor k 2 in our theory 
and —60 times the result of another closest rival theory of 
Spyrou and Papadopoulos [3l| . Therefore, in a way this 
theory is in remarkable agreement with that of Blanchet 
and Schafer. As of today this value remains below the 
accuracy in the measurement of ft [8j, |29j, l3fj . The slight 
difference between the two theories is due to the differ- 
ent methods employed in obtaining the orbital energy 
equations. Eq. (|3.3p in present theory and Eq. (3.35) 
in the theory of Blanchet and Schafer [3j. Present the- 
ory uses basic mathematical tools involving very few ap- 
proximations for deriving the exact orbital energy equa- 
tion, whereas Blanchet and Schafer uses advanced math- 
ematical tools involving plenty of approximations. These 
methods of post-Newtonian wave generation formalism 
are described in Blanchet and Damour [32j , Damour and 
Deruelle [TH, [26|. Again the stress energy tensor used 
by Blanchet and Schafer) uses Newtonian gravitational 
potential and classical kinetic energy where as present 
theory uses relativisitic gravitational potential and rela- 
tivistic kinetic energy. The success of this theory is pri- 
marily due to the introduction of the relativistic mass of 
the orbiting body in the Newtonian theory of gravitation. 
This in turn provides further justification for the use of 
relativistic mass. In PR [l| static WEP is modified to 
dynamic WEP which states that the gravitational mass 
is equal to the relativistic mass. 



1 | l /z 1 (l + e)(l-3e) 
4 ac 2 (l - e) 2 



fJ,l((Tl + <7 2 ) 
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1 5(l + e)(l-3e) /2ttG\ 3 . 

1+ 4 c 2 (l- e) 2 rn l{ m 1+m2 ) 3 

(a 1 + <T 2 )m 1 (m 1 + m 2 )~^ 



6c 2 (l - e) 2 (l + (73/24)e 2 + (37/96)e 4 ) 



'2irG 

VftT 



(3.12) 



4. ORBITAL PHASE OF AN INSPIRALLING 
COMPACT BINARY 

In this section, we will simply relate the orbital energy 
Eq. (|2.7p to the 3.5PN accurate gravitational wave-form 
model discussed by Blanchet et al.[f|. Since the orbits 
of the inspiralling compact binaries are circularized, for 
e = 0, Eq. (|2~7) reduces to 



E° = - 



Vim 2 
2a 



1 Ml 1 
1 + -— H — H 

4 ac 2 8 a 2 c 4 64 a 3 c 6 



f'l 



(4.1) 
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We have the post-Newtonian parameter 7 given by, 



CONCLUSION 



7 



Gm G(nii + TO2) 



/'-2 



r^c ac* ac^ 

Substitution of Eq. P~2]) in Eq. fl~Tl) yields, 



( 4 - 2 ) 



E° = - 



Mc 2 7 



4 Vm 2 / 8 \m 2 ' 



64 V 7Tl2 



1 



(4.3) 



Equation Eq. (|4.3p is comparable to those given in [5|, 
7] . The rest of the procedure for computing the orbital 
phase <p in terms of the frequency related parameter x 
can be same as described in E|. One can either use 
the same 3.5PN accurate expression for the gravitational 
wave flux £ [3, la] or reevaluate C in consideration of the 
new equations of motion resulting from Eq. . 



We have derived precise value for the orbital period 
derivative which is very close to the currently accepted 
theoretical value derived by Blanchet and Schafer 
This needs to be experimentally verified for establishing 
the correct theoretical approach. It would be impossible 
to derive the expression for the orbital period derivative 
in this theory without the use of relativistic mass for the 
orbiting body in the inverse square law of gravitation. In 
PR [l|] static WEP is modified to dynamic WEP which 
states that the gravitational mass is equal to the relativis- 
tic mass. It is the use of relativistic mass that eliminates 
the need for higher order xPN corrections to the orbital 
energy part of the balance equation. Further justifica- 
tion for using the relativistic mass is discussed at length 
in 0. 
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